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1 . REVIEW  OF  ACCOMPLISHMENTS  1976  - 1978 


t. 
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i; 


The  highlights  of  accomplishments  during  the  period  are 
the  construction  and  characterization  of  the  turbulence  chamber 
for  optical  transmission/  eind  its  application  to  the  comparison 
of  optical  transmission  using  a corner-cube,  plane  mirror  and 
direct  transmission.  With  the  objective  of  clarifying  the 
theoretical  models  of  optical  transmission  through  the  turbu- 
lent atmosphere,  we  have  built  a turbulence  chamber  and  con- 
firmed its  isotropic  and  homogeneous  nature  of  turbulence  as 
2 2 

well  as  its  C_  and  C , inner  and  outer  scales  of  turbulence 
T n 

and  L^.  Since  the  strength  of  this  laboratory  generated 
2 -11  -2/3 

turbulence  = 3 x 10  m ' is  about  1,000  times  larger 
than  that  of  a typical  real  atmosphere,  it  enables  us  to 
simulate  the  essential  features  of  optical  transmission  under 
the  well-controlled  condition.  The  turbulence  thus  generated 
is  shown  to  be  locally  isotropic  and  homogeneous  within  the 
region  where  optical  transmission  experiments  were  performed. 
The  Kolmogorov-Obukhov  2/3  power  law  for  the  temperature 
structure  fvinction  was  shown  to  be  valid  within  the  localized 
region  of  6 cm  diameter.  Its  details  are  described  in 
Appendix  1. 

The  power  spectrum  of  temperature  fluctuations  in  the 
turbulence  chamber  closely  followed  theoretically  predicted 
-5/3  power  law  for  inertial  sub-range.  From  the  power  spectrum 
measured,  we  obtained  inner  and  outer  scales  of  turbulence 
(i  and  L , respectively),  using  I.  Taylor's  hypothesis 


(2) 


of  frozen  turbulence:  t = 5 mm  emd  L = 6.5  cm. 

o o 

It  was  possible  with  the  turbulence  chamber  to  stvidy  thv? 
optical  transmission  through  turbulence  for  both  regions  with 
Fresnel-zone  size  smaller  and  larger  than  the  inner  scale  of 
turbulence.  We  have  thus  been  able  to  verify  some  of  the 
theoretical  predictions  of  irradiance  fluctuations  for  both 
regions.  The  power  spectrum  of  irradiance  fluctuations  follow 
-8/3  power  on  frequency  in  close  agreement  with  the  theory. 
Details  are  described  in  the  manuscript  to  be  submitted  to 
Applied  Optics.  (Publication  #2) 

By  increasing  the  optical  path -length  through  the  turbu- 
lence chamber#  we  could  reach  the  so-called  saturation  region 
of  irradiance  fluctuations.  The  coefficient  of  variation 
of  irradiance  fluctuations  becomes  1.6#  which  is  close  to  the 
observed  value  in  the' typical  long  path  propagation  in  real 
atmosphere  under  strong  wind  conditions . These  results  will 
be  described  in  the  manuscript  to  be  svibmitted  to  Radio 
Science.  (Publication  #3) 

Comparison  of  a corner-cube  retroreflector  and  a plane 
mirror  in  irradiance  fluctuation  of  the  transmission  through 
the  turbulence  has  indicated  that  the  path  of  propagation 
for  the  plane  mirror  is-  statistically  more  independent# 
whereas  for  the  corner-cube  the  transmission  paths  are  some- 
what correlated.  The  power  spectruros  of  irradiance  fluctua- 
tions tend  to  be  wider  for  the  plane  mirror  case  than  the 
corner-cube  case.  This  can  be  explained  qualitatively  consi- 


dering  the  Doppler  shift  of  the  scattered  light  by  the  tur- 
bulence. The  reflected  beeun  at  corner-cube  has  coinponents 
opposite  to  that  of  incident  with  respect  to  the  transverse 
air  flow.  Thus,  the  laser  beam  of  frequency  up-shifted  by 
the  Doppler  effect  will  suffer  the  down-shift  on  the  way 
back.  On  the  other  hand,  the  beam  reflected  by  a plane 
mirror  will  produce  doppler  shift  in  the  same  direction  as 
the  incident  beam.  This  is  described  in  Appendix  II. 

Single-ended  irradiance  fluctuations  measurement  using 
a comer-cube  re troref lector  is  convenient  from  a practical 
standpoint;  scientists  at  both  Naval  Ocean  Center  and  National 
Weather  Service  have  approached  us  to  provide  further  infor- 
mation on  this  experiment. 

We  have  performed  preliminary  experiments  of  intensity 
correlation  experiments  using  the  64  diode  array  interfaced 
to  PDP  15/20  data  acquisition  system.  (See  Photograph  on 
Page  6.)  This  may  become  a useful  tool  for  measuring  the 
angle  of  arrival  fluctuations  and  also  for  characterizing 
the  saturation  region  in  which  the  multiple  scatterin  effect 
becomes  significant. 


- ^ • ' 1 
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APPENDIX  I 


AN  ATMOSPHERIC  TURBULENCE  CHAMBER  FOR  OPTICAL  TRANSMISSION 
^ EXPERIMENT:  CHARACTERIZATION  OF  THERMAL  METHOD*+ 

Introduction 

Artificially  generated  atmospheric  turbulences  have  been 

utilized  in  a few  optical  propagation  and  imaging  experiments. 

Most  of  them,  however,  were  performed  without  thoroughly  speci- 
e's 

fying  the  strength  of  turbulence.  We  have  designed  and  built 
a sipole  atmospheric  tu.rbulence  chamber  with  the  intention  of 
performing  optical  wave  propagation  experiments  through  the 
well-defined,  homogeneous  and  isotropic  turbulence.  Since  the 
refractive  index  fluctuations  in  the  atmosphere  are  mainly 
dependent  on  the  temperature  fluctuations,®  we  utilized  the 
heater/blowers  for  generating  the  stronger  refractive  index 
fluctuations  and  made  the  turbulence  more  homogeneous  and  iso- 
tropic by  passing  through  metallic  screens  with  random  foils 
and  2 mm  screens.  In  this  paper,  we  shall  characterize  the 
atmospheric  turbulence  in  the  chamber  by  the  temperature 
structure  constant  which  is  determined  from  temperature 

fluctuation  measurements , and  then  derive  the  refractive  index 

structvire  constant  C ^ from  C,.^.  The  of  the  turbulence 

n I I 

chamber  was  found  to  be  1,000  times  stronger  than  that  of 
the  typical  atmosphere.  The  characterization  by  the  irradi- 

ance  fluctuation  measurements  will  be  described  in  the  second 
paper. 

The  histogram,  variance  and  power  spectrum  of  fluctuations 
of  the  temperature  difference  betV'feen  two  micro- thermocouple 
probes  were  obtained.  The  experimental  data  was  processed  using 
a PDF  15/20  system.  The  temperature  structure  function  (the 
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mean  square  temperature  difference  between  two  probes)  was 
measured  as  a function  of  the  probe  separation . It  is  shown 
that  the  temperature  structure  function  depends  on  the  2/3 
power  of  the  probe  seperation  from  2 to  6 cm  in  accordance  with 
the  Kolmogorov-Obukhov  similarity  law.  The  temperature  struc- 
ture constants  were  obtained  from  these  measurements  of 
temperature  structure  function.® 

The  power  spectrum  of  temperature  fluctuations  measured 
inside  the  turbulence ’ chamber  is  dependent  on  the  -5/3  power 
of  the  frequency  within  the  frequency  range  from  5 to  80  Hz. 

This  agrees  with  the  well-known  f“5/3  fQj,  inertial 
S’Jib-range  of  homogeneous,  isotropic  turbulence.®'^®  By 
assuming  that  turbulence  eddies  were  "frozen"  while  flowing 
at  the  mean  velocity  0.41  m/sec,  we  estimated  the  inner  and 
outer  scales  of  turbulence  of  the  upper  and  lower  frequency 
limits  mentioned  above;  Jlo=5  mm  and  Lo=8.2  cm,  respectively. 

The  refractive  index  structure  constants  C..  were  derived 

n 

from  the  temperature  structure  constants  by  using  the  tempera- 
ture dependence  of  the  refractive  index  of  the  atmosphere.  The 
of  our  turbulence  chamber  is  1,000  times  larger  than  that 
of  typical  turbulence  atmosphere. 

Different  methods  for  generating  the  homogeneous  isotropic 
turbulence  are  compared  to  our  method.  Possible  improvements 
of  our  turbulence  chamber  and  some  meaningful  aerodynamic  1 
experiments  are  discussed. 

iiic  measurement  or  irequency  response  function  of  micro- 
thermocouples is  described  in  Appendix  A and  the  rms  temperature 
gradient  and  characteristic  temperatxire  of  the  turbulent 
temperature  field  were  s aluated  by  using  experimental  results 
and  is  described  in  Appendix  B. 


The  design  of  the  turbulence  cheunber  is  schematically 
illustrated  in  Fig.  1.  The  turbulence  chamber  with  internal 
dimensions/  2.69  m long/  0.78  m wide  and  0.23  m high/  was 
assembled  on  a shock-mounted  optical  bench.  The  internal 
surface  of  side-walls / top  covers  and  aluminum  plates  of  the 
op.tical  bench  were  coated  by  a special  black  paint  (3M  101 
CID)  in  order  to  minimize  disturbances  due  to  the  reflection 
of  the  scattered  light  by  the  chamber  walls. 

Ten  small  electric  heater/blowers  (Sears / 764/  R,  Model 
344)  were  installed  on  the  one  side  of  the  turbulence  chamber 
The  turbulent  flows  . generated  by  the  heater/blowers  are  . ; 

parallel  to  the  optical  bench  surface  and  pe’.^andicular  to 
the  optical  transmission  direction.  The  height  of  the  center 
of  the  flow  was  almost  equal  to  our  .standard  height  for 
optical  measurements/  0.11  m.  The  turbulent  flow  generated 
by  the  heater/blowers  was  homogenized  by  passing  through  a 
screen  of  aluminum  foil  (commonly  used  in  electronic  instru- 
ments for  ventilation  filters)  and  then  through  three  layers 
of  2 mm  square  aluminum  meshes,  each  of  which  is  separated 
by  1 cm.  The  turbulent  fl6w  produced  by  heater/blowerr 

s. 

operated  at  the  rate  of  750  watts  was  found  convenient  for 
the  various  optical  transmission  experiments , because  of  its 
large  strength  of  turbulence. 


(9) 

Average  Wind  Velocity  and  Temperature  Distribution 

The  wind  velocity  and  temperature  were  measured  along  the 
flow  direction  (x)  by  a hot  wire  anemon-r-ter  (Alnor  Thermo- 
Anemometer  8500)  throughout  the  turbulence  chamber  at  t’ie 
standard  height  0.11  m.  The  results  of  wind  velocity  measure- 
ments are  shown  in  Fig.  2.  The  uniformic/  of  the  wind  velocity 
in  the  flow  direction  (x)  improved  with  distance  from  the 
last  screen.  At  x=*0.61  m,  the  wind  velocity  became  fully 
uniform.  The  average  wind  velocity  in  the  r2Uige  0.30  m to  0.61 
m was  0.41  m/sec  at  the  standard  height  of  0.1103  m. 

The  temperatxire  profile  for  the  case  in  which  the  heater/ 
blower  was  set  to  750  watts  is  shown  in  Fig.  3.  The  temperature 
averaged  over  the  chamber  was  53“  C.  Comparing  the  results 
in  Figs.  2 and  3,  we  notice  that  the  te'mperature  distribution 
became  viniform  at  a much  faster  rate  than  the  wind  velocity 
distribution. 

The  development  of  turbulent  flow  inside  the  chamber  may 
be  explained  in  the  following  way . ^ ^ A flow  generated  by  each 
heater/blower  at  the  inlet  of  the  chanber  forms  boundary  layers 
on  both  upper  and  lower  walls,  owing  to  the  viscous  friction. 

As  flow  propagates  along  the  x direction,  the  interaction 
between  flow  shear  and  inlet  turbulence  increases.  Thus,  the 
thickness  of  boundary  layers  formed  on  both  walls  will  increase 
in  the  downstream.  As  a result,  the  uniformity  of  the  bulk 
flow  in  the  y direction  will  improve.  This  is  exhibited  by 
Fig.  4,  where  wind  velocity  profiles  along  the  y axis  (height) 
for  X = 0.013  m,  0.30  m and  0.61  m from  the  last  screen  at  the 
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distance  z»1.09  m are  illustrated.  It  was  concluded  that  a 
region  in  the  chamber  of  about  0,31  m wide  (x) , 0.076  m high 
(y)  and  2.54  m long  (z)  has  relatively  uniform  velocity  and 
ten^erature . 

The  Reynolds  number  Re=ud/v  for  the  region  of  uniform 
flow  in  the  turbulence  chamber  is  Re=6287,  where  the  average 
flow  velocity  u=0.41  m/sec,  distance  between  walls  d=0.23  ro 
and  the  kinematic  viscosity  of  air,  v=1.5  x 10” 5 m/sec^  are 
used..  According  to  the  standeurd  reference  books  by 
Schlichtihg-  and  Vennard  and  Street, the  critical  Reynolds 
niambers  for  flow  between  parallel  walls  are  2740  and  1000, 
respectively.  Since  the  Reynolds  nvimber  6287  calculated  for 
our  turbulence  chamber  is  much  greater  than  these  critical 
numbers,  the  turbulent  flow  within  this  region  of  the  chamber 
should  be  locally  homogeneous  and  isotropic. This  is  well- 
evidenced  by  the  experimental  results  described  in  the  next 


section. 


1 
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Measiareiaents  of  Temperature  Structiure  Functions 

The  ten^jerature  structure  ftanction  Py(r)  is  defined®  as 

the  mean  square  of  the  temperature  difference  between  tvio 

points,  r and  r : 

1 2 

D^-Cr  , r ) = < ET(r  ) - T(r  ) ]2>  . (1) 

1 1 2 -1  2 

When  the  temperature  field  is  homogeneous  and  isotropic,  the 

stmicture  fxiuction  depends  only  on  the  distance  between  two 

points  r = I r - r I : 

12 

(r  , r ) = Dt-  (r)  ,(2) 

The  experimental  set-up  for.  the  temperature  structure 
function  measurement  is  illustrated  in  Pig.  5.  The  disteuice  r 
between  two  micro-thermocouple  probes  mounted  on  an  electro- 
mechanical translating  unit  can  be  adjusted’ from’ 2 mm  up  to 
20  cm  by  the  increment  of  10  jim.  The  mechanical  translation 
unit  was  driven  by  an  electrical  stepping  motor  (SLO-SYN 
SS50-1009,  Superior  Electric  Co.)  with  an  external  controller, 
which  can  be  triggered  by  the  pulse  from  the  PDP  15/20  data 
acquisition  system  through  a digital-to-analog  convertor. 

This  electro-mechanical  translation  unit  was  placed  on  the 
top  cover  of  the  turbulence  chamber.  Only  the  thesnnocouple 
probes  supported  by  the  vertical  stainless  steel  tubes  were 
exposed  to  the  turbulent  flow  (See  Fig.  1) . The  thermocouple 
probes  project  horizontally  up-stream  by  8.5  cm  from  the 
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vertical  supporting  tubes.  The  temperature  fluctuations  in 
the  turbulent  air  flow  were  measured  by  using  differential 
micro-thermocouples  of  ohromel-constemtan  (Type  E,  Hy-Cal 
Engineering)  , of  which  the  di.'ameter  is  7.6  ym  emd  the  speci- 
fied time  constant  is  less  than  0.1  m/sec.  The  temperature 
voltage  conversion  factor  K of  the  micro-thermocouples  was 
calibrated  by  using  two  electric  ovens  while  keeping  the 
temperature  of  one  junction  at  53*  C and  varying  the  tempera- 
ture of  the  other  oven:  K = 0.07184  x 10“ ^ volt/®K. 

The  frequency  response  of  micro-thermocouples  is  constant 
from  DC  to  200  Hz  and  then  decrezuses  to  % of  the  DC  level 
at  770  Hz  (See  Appendix  A) . Therefore,  the  micro- thermocouple 
system  can  accurately  measure  the  temperature  fluctuations  of 
the  turbulent  atmosphere,  of  which  the  predominemt  frequency 
component  is  limited  within  200  Hz. 

The  output  voltages  of  two  thermocouples  were  eunplified 
by  differential  amplifiers  (Tektronix  AM  502),  respectively, 
and  two  single-ended  outputs  from  these  amplifiers  were 
connected  to  the  differential  input  connectors  of  the  third 
Tektronix  amplifier  (AM  502)  through  two  50  n coaxial  cables 
of  equal  length.  The  temperature  difference  at  two  points 
can  be  derived  from  the  differential  voltage  Uo  measured  by 
using  the  formula 

T -T  = ^ u„  , (3) 

^ ^ G G G 

1 2 3 

where  K the  thermocouple  constant,  G , G and  G are  the 

12  3 


gains  of  three  amplifiers  mentioned  etbove,  u«  is  the  voltage 
measured  by  an  analog- to-digital  convertor.  By  inserting 
Eq.  3 into  Bq.  1,  we  obtadn  the  temperature  structure  function 
expressed  in  terms  of  quantities  measured  and  calibrated  as 
follows : 

<(T  - T )^>  = ^ <o,2>  , (4) 

1 2 (G  G G )2 

12  3 

where  the  veuriance  <Uo^>  of  the  output  voltage  from  the  third 
eunplifier  was  determined  from  the  digitized  output  voltages  ‘ 
by  the  statistical  data  acquisition  system  illustrated  in 
Pig.  6. 

In  the  so-called  inertial  sxab-range  of  the  hompgeneous 
isotropic  turbulence,  the  teir^erature  structure  function  will 
depend  on  the  2/3  power  of  distamce  r between  two  probes:® 

Dy  (r)  = j-2/3  ^ (5) 

where  the  distance  between  probes  should  be  much  larger  than 
the  inner  scale  of  turbulence  lo  but  much  smaller  than  the 
outer  scale  of  turbulence  L®:  <<  r<<  Lp/  and  is  the 

temperature  structure  constant.  This  2/3  power  equation  is 
often  called  the  Kolmogorov-Obukhov  similarity  law. 

The  temperature  structure  functions  at  a location  x = 0.38  m 
y = 0.12  m and  z = 1.09  m,  in  the  turbulence  chamber,  measured 
as  a function  of  the  probe  seperation  r,  are  illustrated  in  the 
logarithmic  scales  in  Pig.  7.  By  inspecting  these  results, 
we  notice  that  for  the  probe  seperation  r between  2 and  6 cm, 
the  2/3  law  holds  reasonably  well.  The  data  points  for  the 


distance  r ssaaller  th^  2 cm  deviate  from  the  2/3  line.  How** 
ever,  the  experimentai  results  do  not  agree  with  the  theoretical 
ten?)erature  structure  function  for  distance  within  the  inner 
scale  of  t\irbulence:  0^  (r)  «r^.  This  discrepancy  is  most 
likely  due  to  the  disturbance  produced  by  stems  holding  a 
micro-thermocouple . The  structure  function  measured  starts 
to  deviate  from  the  2/3  law  at  the  probe  distance  r = 6 cm 
and  attains  a constant  value  for  r = 7 cm  and  higher.  This 
ceui  be  explained  by  considering  that  for  a large  probe  sepera- 
tion  larger  than  the  outer  scale  of  turbulence,  the  tempera- 
ture fluctuations  at  the  two  probes  are  uncorrelated.  Since 
the  cross  term  of  temperatures  in  expansion  of  <(T^  - 
veuiish,  temperature  structiire  function,  D^(r)  becomes  inde- 
pendent of  r.  Furthermore,  the  measured  Dy(r)  indicates 
that  the  outer  scale  of  turbulence  is  of  the  order  of  7 cm 
and  the  temperature  field  in  this  turbulence  chamber  is 
locally  homogeneous  and  isotropic  within  the  6 cm  diameter 
at  the  standard  height. 

Using  the  log  D^(r)  extrapolated  by  the  2/3  line  at 
r = 1 in  Pig.  7,  we  obtain  the  strength  of  temperature  fluctua- 
tions = 52.9  m“V3,  Frcm  the  measured  we  also 

O 

derived  the  rms  temperature  gradient  1.898  K/cm  and  the 

O 

characteristic  temperature  0.949  k {See  appendix  B) . 

The  uniformity  of  the  strength  of  turbulence  within  the 
chamber  was  evaluated  by  performing  a series  of  measurements 
of  mean  square  temperature  difference  at  a probe  seperation 
2.89  cm  where  the  2/3  law  approximately  holds. 


From  the 
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measured  veuriances  of  temperature  difference  fluctuations,  we 
obtadned  the  defined  above.  Ihe  measured  are  shovm 
in  Fig  8.  It  should  be  noted  that  has  a maximum  value 
down  the  turbulent  flow  from  the  last  screen  while  both  the 
average  velocil^  and  temperature  distributions  will  monotomi- 
cally  decrease.  This  implies  tdiat  tlie  mean  square  temperature 
gradient  takes  the  maximum  value  somewhat  down  the  flow  around 
X = 0.46  m (refer  to  Eqs.  (A'3)  and  (A4)  . 

The  turf)ulence  strength  parameter  of  temperature  fluctua- 
tions is  small  neeir  the  inlet  of  the  chamber  (i.e.,  just 
near  the  screens  adjacent  to  the  heatsr/b lowers.)  This  indi- 
cates that  the  ten5)erature  fluctuations  of  grid-generated 
turbulence  is  small.  As  the  heat  flows  along  x direction, 
the  small  eddies  interact  with  eddies  generated  in  the  shear 
layer  adjacent  to  constant  temperature  walls.  As  eddies  from 
both  sources  mix,  tenperature  fluctuations  increase.  In  the 
region  between  x = 0.30  m and  x = 0,61  m,  the  interaction 
between  the  two  sources  of  eddies  becomes  the  largest.  In 
the  further  downstream,  the  effects  of  wall  temperature  become 
dominant  and  the  net  tenperature  fluctuations  decrease. 
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Power  Spectnaa  of  Temperature  Fluctuations 

By  calculating  the  time  average  of  absolute  value  squued 
of  Fourier  trams form  of  digitized  thermocouple  output  voltage / 
we  estimated  the  power  spectrum  of  temperature  fluctuations. 

512  digital  sanples  with  sauipling  interval  1 m sec  were 
processed  by  the  Fast  Fourier  Transfor^n  (FFT)  method  and  were 
averaged  over  20  results.  The  power  .spectrum  of  temperatiire 
fluctuations  at  the  location  x =»  0.38,  y » 0.11  n,  z * 1.09  m 
and  the  probe  distamce  r = 2.89  cm  acre  illustrated  in  Fig.  9. 
Over  the  frequency  region  from  f = 10  Hz  to  30  Hz,  the  esti- 
mated power  spectrum  W(f)  is  proportional  to  f**®/^.  At  the 
frequency  higher  than  30  Hz,  the  power  spectrum  approximately 
follows  the  f”^  law,  which  was  obtadned  by  W.  Heisenberg®  in 
the  viscous  dissipation  region.  At  the  frequency  region  below 
10  Hz,  it  approaches  a constant  value  in  a manner  similar  to 
the  atmospheric  turbulence  model  disscussed  by  Strohbehn. 

In  order  to  determine  the  transition  frequency  from  the 
inertial  sub-range  to  the  viscous  dissipation  range,  we  also 
calculated  f^  W(f)  and  showed  in  Fig.  9.  The  upper  and  lower 
crossing  points  of  the  f”®/®  line  for  inertial  sub-range  with 
the  f"^  and  constant  lines  in  Pig.  9 are  given  by  fj^  ==  80  Hz 
and  fj^  = 5 Hz,  respectively.  Since  the  inner  scale  of  txirbu- 
lence  equal  to  Ao  of  the  inertial  sub-range  is  given  by)lo  = u/f^, 
by  inserting  the  average  wind  velocity  u = 0.406  m sec  and  f^^  = 
80  Hz,  we  obtain  Ao,  = 5 mm.  Likewise,  by  inserting  the  same 
average  wind  velocity  and  f = 5 Hz  into  the  outer  scale  of 
turbulence  L,  =u/f.,  we  obtain  the  outer  scale  of  turbulence 

Xf 

Lo  = 8 cm.  The  Reynolds  numbers  for  these  scales  of  turbulence 
Re^  = uio/v  and  Re^  = uL„/v  are  given  by  Re*’  = 135  and  Re^  = 2, 
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165,  respectively. 
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Statistics  cf  Tenperature  Fluctuationg 

The  512  channel  histogram  of  tenperature  difference  at 
r » 2 cm  were  obtained  by  using  the  POP  15/20  data  acquisition 
system  for  vzurious  locations  of  the  turi>ulence  chamber.  From 
these  histograms  r we  derived  the  central  moments  p.  amd  cumu-* 

m 

lamts  Xj,  up  to  the  eighth  order. The  coefficients  of 

variation  y / skewness  y and  excess  y are  calculated  from 
• 12 

measured  moments  and  cximulamts,  including  their  error 
estimates : 


Y * <t/u  , (6) 

« 1 

Y » p /a^,  (7) 

1 3 


amd 

yap  /a'»  - 3 , 

2 4 


(8) 


where  a is  the  standaurd  deviation  and  is  the  variance.  The 
smoothed  frequency  distribution  curves  of  tenperature  can  be 
obtained  by  inserting  the  above  coefficients  into  the  Edgeworth 
series : ^ ^ 

f (5)=  ♦(?)  - Yj  (5)  + yj  (5)  + (9) 


Y^2  ^(6)  (5)  + ... 

where 

1 

^ (g)  = — exp  (-  ■5'  ) f 

/'  2ir 

and 

ia  the  n-th  derivative  of 

(5)  = (-1)”  4.(5) 


f 


(10) 
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(x-Pi)/a  . 


(i: 


(5)  is  the  Hermite  polynomial  of  degree  n: 

H,  (x)  » 1,  Hi  (x)  * X,  H2  (x)  * x^  - 1 

H3  (x)  * x^  - 3x,  H4  (x)  » x**  - 6x^  + 3 . (i; 

The  smoothed  frequency  distribution  curves  for  a few  locations 
along  the  air  flow  are  shown  in  Pig.  10.  The  coefficients  of 
variation,  skewness  and  excess  are  shown  as  a function  of  x 
in  Pig.  11. 

According  to  these  results,  the  frequency  distribution 
of  ten^erature  fluctuations  near  the  screens  is  very  close  to 
the  gaussian  distribution.  {Senerally,  in  the  gaussian  random 
process,  the  different  frequency  spectral  components  are  sta- 
tistically independent.  However,  when  the  different  frequency 
components  eure  not  independent  of  each  other,  statistical 
distributions  become  non-gaussian.  In  the  turbulent  flow  of 
inertial  sub-range,  lower  frequency  spectral  components  will 
transfer  energy  into  the  higher  frequency  components.  Hence, 
the  histograms  of  temperature  fluctuations  in  the  turbulent 
flow  are  generally  non-gaussi&n.  In  the  region  of  maximum 
coefficient  of  variation,  x=  0.46  m,  the  wide  range  of  frequency 
spectral  components  representing  various  sizes  of  eddies  should 
be  interacting  most  actively.  According  to  Pig.  11,  the  excess 
coefficients  increase  along  the  flow  and  becomes  maximum  at 
x=  0.54  m.  However,  the  coefficients  of  skewness  increase 
raonotomically  with  distance. 
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Derivation  of  C.  from  C,. 

' ' ■ n ""  I 

The  refractive  index  fluctuations  of  the  homogeneously 
ttirbulent  atmosphere  at  optical  wavelength  are  mainly  due  to 
the  temperature  fluctuations.  The  refractive  index  structure 
function  D„(r)  » <[n(r  ) - n(r  )]^>  in  the  inertial  svib-range 
cam  be  expressed  in  the  same  2/3  law  as  the  temperature 
structure  function: 


D_(r)  * C 2 r2/3. 
n n 


(14) 


The  refractive  index  structure  constant  can  be  expressed 
in  terms  of  the  temperature  structure  constant  by  using  the 
equation  of  refractive  index  as  a function  of  ten^erature  and 
pressure  with  wavelengtli  as  a parameter.  The  refractive  index 
of  the  air  at  the  temperature  T*K,  the  pressure  p mb  and  the 
wavelength  X ym  can  be  approximated  by  the  following  formula: 


77.6  P 
T 


0,^0753  1 ^ 10-6  . 

X2  J 


(15) 


Since  the  refractive  index  under  the  constant  pressure 
depends  only  on  the  tenperature,  we  obtain:^® 

C = 77.6  P r ^ ^ 0.00753 

n t2  L X2 

For  the  normal  atmospheric  pressure  p * 1,013.25  mb  and  the 

wavelength  X = 0.6328  wm  and  the  average  temperature  of  the 

turbulent  atmosphere  T = 326.16  “K,  we  obtain  = 0.75302  x 

10"®  C,  and  C 2 = 0.56704  x 10" ^2  q 2 ^ 

I n I 
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Discussion 

In  the  abo've  derivation,  we  assumed  only  the  ten5>erature 
dependence  of  refractive  index  of  the  atmosphere  and  neglected 
its  pressure  dependence.  However,  we  may  have  to  include  some 
effect  of  pressure  fluctuations.  Since  in  case  of  sound  waves, 
pressure  effect  is  usually  treated  as  an  adiabatic  process,  it 
would  be  rezisonable  to  assume  that  the  pressure  fluctuations 
induce  the  temperature  fluctuations  through  the  adiabatic 
thermodynamical  process.  Bien,  asstaming  the  refractive  index 
n is  a function  of  the  temperature  T and  pressure  P and  using 
a relation  for  the  adiabatic  process 


^ _ _x__  2, 

dT  Y~1  T 


Cp/Cy)  (17) 


we  obtain 

(adiabatic) 

n 


77.6  P 

t2 


_1_  'i  + 1v°°753  x10“&  C, 

Y-1  I 


(18) 


Inserting  y = 1.40  of  the  normal  air,  we  obtain  (adiabatic) 

= 2.48  C . Thus,  it  is  essential  for  establishing  the  atmos- 
pheric  model  to  compare  the  directly  measured  by  the  optical 
method  and  the  obtained  from  Ct  . 

When  the  turbulence  of  the  atmosphere  is  neither  isothermal 
nor  adiabatic  process,  we  may  have  to  include  the  pressure 
structxore  constant  Cp  and  the  pressure  and  temperature  cross- 
structure  C - . In  case  of  the  atmospheric  turbulence  in  the 
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outside  field,  we  mu^t  also  include  the  effect  of  the  humidity 
fluctuations  and  the  density  of  other  gases  such  as 

We  would  like  to  point  out  two  interesting  research  areas 
pertinent  to  our  turbulence  experlir^nt,  although  they  are  not 
directly  related  to  the  optical  propagation.  First,  we  may  be 
able  to  clarify  the  physical  meaning  of  non-gaussian  sta;istics 
of  temperature  fluctuations  exemplified  in  Figs.  10  and  11,  by 
using  some  of  the  recent  research^  results-^*’  on  the  tenperature 
fluctuations.  Second,  by  measuring  the  power  spectrum  of 
temperature  fluctuations  within  the  turbulent  boundary  layer, 
we  may  be  able  to  determine  the  model  of  turbulence  in  the 
viscous-convective  range.  This  is  an  interesting  subject, 
because  the  "bump"  in  the  power  spectrum  of  temperature 
fluctuation,  that  is,  the  deviation  from  the  -5/3  power  law, 
have  been  observed  in  the  atmospheric  surface  layer^^  and  is 
explained  by  the  -1  power  law  for  the  viscous-convective 
range .22,23 

Two  other  methods  for  generating  ^e  isotropic  and 
homogeneous  txirbulence  have  been  reported;  a "porcupine"  like 
array  of  jet  streams  used  by  Betchov  and  Lorenzen^**  and  the 
heated  grid  in  a supersonic  wind  tunnel  used  by  Shih-chun  Lin 
and  Shao-chi  Lin. 25 '26  compared  to  these  methods,  our  turbu- 
lence chamber  is  much  simpler  to  construct  and  more  convenient 
to  apply  to  the  optical  transmission  measurements. 

The  performance  of  our  turbulence  chamber  may  be  improved 
by  some  modifications;  for  instance,  by  adding  the  array  of 
exhaust  fans  at  the  outlet,  we  should.be  able  to  widen  the 


- ■ 

region  of  uniform  strength  of  turbulence.  It  is  also  desirable 
to  channel  the  outlet  flows  to  the  air  exhaust  system  in  the 
laboratory,  because  even  higher  tenperature  setting  of  heater/ 
blowers  may  be  utilized.  By  making  an  array  of  hol?s  on  the 
walls,  we  may  enhance  tenperature  mixing  and  generate  the 
stronger  turbulence. 

After  the. manuscript  of  this  paper  was  prepared,  we 
noticed  an  interesting,  recent  paper  by  Bissonnette^^  on  the 
optical  transmission  in  the  aqueous  turbulence  with  the 
specified  strength  of  turbulence,  ihe  objective  of  the 
laboratory  simulation  described  in  the  paper  is  very  similar 
to  ours.  Our  atmospheric  turbulence  chamber,  however,  has 
advantages  such  as  convenience  of  inserting  various  optical 
and  aerodynamical  instmaments  in  the  chamber  coiipared  to  the 
water  chamber  used  by  Bissonnette. 

Conclusion 

The  average-  strength  of  refractive  index  fluctuations 
of  the  turbulent  flow  generated  within  the  chamber  was  10 ^ 
times  larger  than  that  of  the  typical  atmosphere  in  the  field. 
The  turbulence  was  locally  homogeneous  and  isotropic  within 
the  region  of  6 cm  diameter.  The  average  strength  of  turbu- 
lence  was  adequately  homogeneous  for  the  purpose  of  optical 
transmission  measurements  within  the  region  0.5  m x 0.05  m x 
2 m at  the  standard  optical  height.  Characteristic  features 
of  the  generated  turbulence  are  reproducible  within  the 
experimental  errors . Various  theoretical  models  of  optical 
wave  propagation  in  the  turbulent  atmosphere  are  being  tested 
by  using  the  laboratory  generated  turbulence  with  defined 
C 2.28  determination  of  C„^  by  irradiation  fluctuation 

measurements  will  be  described  in  a seperate  paper. 
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Appendix  A:  Meaisurement  of  Frequency  Response  of  a Micro- 
Thermocouple  by  Using  a Modulated  Laser  Beam. 

A modulated  50  mW  He-Ne  632.8  nm  laser  beam  (single-mode 
Spectra-Physics  125  A with  Fedjry-Perot  Etalon  in  the  cavity 
extender  and  Spectra-Physics  320  Modulator)  was  focused  on 
micro-thermocouple  junctions  while  the  other  junction  was 
covered  by  a plastic  case.  The  frequency  response  of  the 
differential  micro-thermocouple  system  was  measured  by  t2dcing 
the  ratio  of  the  output  signal  from  the  thermocouple  system 
with  those  from  a photomultiplier  tube  which  detects  the 
modulated  laser  beam  directly.  In  order  to  determine  the 
frequency  response  xinder  the  operating  condition,  the  measure 
ments  were  perfoimied  with  heated  air  flow.  The  measured 
frequency  response  is  shown  in  Pig.  Al. 
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J^pendix  B:  Evaluation  of  Average  Temperature  Gradient  emd 
Characteristic  Temperature  From 

The  temperature  structure  constant  Cj^  can  be  expressed  as® 

, „2^/e^l/3  (Al) 

where  rhe  parameter  a is  given  by  the  average  energy  dissi- 
pation rate  e and  the  kinematic  viscosity  v as: 

a-  ej/3v  . (A2) 

A parameter  N,  which  characterizes  the  intensity  of  the 
temperature  fluctuations,  is  given  by  the  thermal  diffusivity 
X and  tlxe  mean  square  temperature  gradient  as 

N =X<  { 7T.)  2>  , (A3) 

and 


X=  k/pO^  (A4) 

P 

where  < is  the  thermal  conductivity,  p is  the  density  of  the 
atmosphere  and  is  the  constant  pressure  specific  heat. 

By  inserting  Eqs.  {A2)  and  (A3)  into  Eq.  (Al) , we  obtain 

Cj^  = X®(AT)^  / (9v^  • 


(A5) 
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The  average  thermal  energy  die<slpation  rate  is  given  by 

ej=  xVa*"*.  (A6) 

Since  we  determined  Cj^  •-  52'.9*C/m'*/^  and  A,  = 0.5  cm, 
inserting  x = 0.200  cmVsec  and  v = 0.1511  cmVsec,  we  obtain 
rms  temperature  gradient 

/ <(7T)  ^>  * 1.898  Vera.  (A7) 


G 


G 


The  thermal  diffusivity  x was  obtained  by  inserting  ic  = 2.6  x 
lO"**  Watt/(cm® ,c) , p = 1.2928/  g/liter  and  * 0.2398  Cal/ 
(g.«®c)  into  Eq.  (A4) . By  using  Eq.  (A6)  and  the  thermal 
diffusivity  x*/  we  obtain  average  thermal  energy  dissipation 
rate  e^: 

Cj  = 0.128  cmVsec^  . (A8) 


The  characteristic  tenperature  of  the  turbulent  tempera- 
ture field  is  defined  by  the  rms  tenperature  gradient  multiplied 
by  the  inner  scale  of  turbulence  Ig  (Ref.  6,  p.  64).  Using 

' ■?' 

5 Ao  = 5 mm,  we  obtain  characteristic  temperature  of  our  turbu- 

E lence  chamber  T = 0.949  ®K. 
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Fig.  3.  Average  air  teiq)erature  within  the  turbulence  chamber  la  shown  as  a function 
of  the  distance  x from  the  last  screen  along  the  direction  of  air  flow. 


(34) 


Schematic  diagram  of  measurement  system  of  tenq>erature 
structure  function  by  micro- thermocouples . 


PC'WER  SPECTRA  OF  TEMPERATURE  FLUCTUATION 
W(f)  and  f^W(f)  (ARBITRARY  SCALE) 


Fig.  9.  Power  spectrum  of  temperature  fluctuation 

2 

W(f)  and  f W{f)  are  shown  as  a function  of 
frequency  (Hz)  where  slopes  of  respective 
regions  are  indicated. 
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Fig.  11.  Coefficients  of  variation,  skewness  and  excess  of  temperature  fluctuations 

are  shown  as  functions  of  the  distance  x from  the  last  screen  along  the 
direction  of  turbulent  air  flow. 
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APPENDIX  II 

i 

COMPARISON  OF  A CORNER-CUBE  REFLECTOR  AND  A PLANE  MIRROR 
. ' 

IN  FOLDED  PATH  AND  DIRECT  TRANSMISSION 
THROUGH  ATMOSPHERIC  TURBULENCE 

INTRODUCTION 

The  objective  of  the  experiment  was  to  clarify  the  different 
characteristics  of  folded  transmission  by  a corner-cube  reflector 
and  a plane  mirror  through  the  atmospheric  turbulence/^  and  also 
to  establish  theoretical  models  for  characterizing  the  irradiance 
fluctuations  measured  for  these  cases.  Due  to  the  unpredictable 
conditions  of  the  real  atmosphere,  we  have  emphasized  on  the  ex- 
periments performed  by  using  an  atmospheric  turbulence  chamber, 
which  can  generate  a locally  isotropic  and  homogeneous  turbulence  | 

I 

in  the  laboratory.  The  results  of  the  field  experiments  is  also  ' 

described. 

) 

EXPERIMENTAL  SET-UP  FOR  THE  LABORATORY  EXPERIMENTS 

The  experimental  set-up  is  shown  in  Fig.  1.  A collimated 

i 

I 

laser  beam  with  the  diameter  of  cross-section  2.54  cm  from  a j 

single-mode  He-Ne  gas  laser  (Spectra  Physics  13.9,  6328  A)  was  j 

1 

r 

used  for  the  laboratory  experiments  using  an  atmospheric  turbu-  1 

! 

lence  chamber.  A corner-cube  reflector  with  the  diameter  , 

2.22  cm  and  a front  surface,  aluminized  plane  mirror  of  the  size 
4.45  cm  X 2.54  cm  with  X/4  flatness  were  used  for  these  experi- 


ments . 


The  light  beams  transmitted  through  the  turbulent  atmos- 
phere were  detected  by  a photomultiplier  tube  (IP21)  with  an 
entrance  diaphragmof  one  mm  diameter,  and  the  output  voltage 
from  the  detector  was  amplified  by  two  cascaded  DC  an?)lifiers 
and  then  was  converted  into  the  digital  form  by  an  Analog-to- 
Digital  convertor  (Preston  Scientific,  14  bits,  plus  sign  and 
100  KC  sampling  rate.)  The  digital  signals  thus  obtained  were 
processed  by  a digital  data  acquisition  system  (PDP  15/20  with 
24  K core  memory.)  The  temporal  frequency  spectra  ox  irradi- 
ance  fluctuations  were  obtained  by  using  the  Fast  Fourier 
Transform  method  and  histograms  of  irradiance  fluctuations 
were  obtained  and  various  moments,  central  moments  and  cumu- 

lants  were  determined  from  these  histograms  (See  Fig.  2.)  The 

2 

turbulence  chamber  with  internal  dimensions,  2.59  m long, 

0.78  m wide  and  0.23  m high  was  assembled  on  a shock-mounted 
optical  bench  (See  Fig.  3.)  The  heated  turbulent  flow  was 
generated  bi  ten  small  electric  heaters  with  blowers.  The 
average  transverse  wind  velocity  at  the  standard  optical  height 

O 

0.12  m was  0.41  m/sec,  and  the  average  temperature  was  53  C. 

The  Reynolds  number  of  this  chamber  was  estimated  as  6287 

with  the  separation  between  upper  and  lower  plates  0.23  m and 

-5  2 

the  kinematic  viscosity  of  the  air  v = 1.5  x 10  m /sec.  The 
value  exceeds  the  critical  Reynolds  numb*  rs  for  the  flow  be- 
tween parallel  walls.  The  turbulence  geiu  rated  was  locally 
homogeneous  and  isotropic  within  a sphere  of  6 cm  diameter.  The 
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average  strength  of  turbulence  was  adequately  homogeneous  for 

this  optical  measurement  within  the  region  0.5  x 0.005  x 2 m 

at  the  standard  optical  height.  The  average  strength  of  re- 

2 -11  -2/3 

fractive  index  fluctuations  was  2.65  x 10  m . The 

inner  scale  I and  outer  scale  of  turbulence  L were  estimated 
o o 

as  1 = 5 mm  and  L = 8 cm,  respectively, 

o o 

RESULTS  OF  LABORATORY  EXPERIMENTS 

The  normalized  power  spectrum  of  irradiance  fluctuations 
measured  was  flatter  and  broadened  for  the  corner -cube  re- 
flector and  was  sharper  for  the  plane  mirror  reflector.  The 
power  spectrum  of  irradiance  fluctuations  in  the  plane  mirror 
reflector  was  broader  than  that  of  the  corner-cube  reflector 
case.  The  histograms  of  corner-cube  reflector  are  closer  to 
the  direct  transmission  of  the  same  transmission  path  length. 
The  statistical  quantities  measured,  the  coefficients  of  varia- 
tion (y  ) , skewness  (y,)  as  well  as  the  quantity  f /2r  XlT  / V 
and  frequency  bandwidth  Af  are  summarized  in  Table  1,  where 
the  critical  frequency  f^  = V / / 2-n  XL,  V the  average  trans- 
verse wind  velocity,  X the  wavelength  of  light  and  L the  path 
length.  The  histograms  obtained  are  shown  in  Fig.  4 

FIELD  EXPERIMENTS 

A similar  experiment  was  carried  out  in  a real  turbulent 
atmosphere  (one  way  path  length  274  m)  on  the  night  of  Dec.  14, 
1976.  The  wind  velocity  was  approximately  3 miles  per  hour. 
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The  results  are  siuimtarized  in  Table  2. 

The  profile  of  the  height  of  transmitted  beam  is  illus- 
trated in  Fig.  5. 

DISCUSSIONS 

The  experimental  results  of  histogreims  may  be  qualitively 
explained  by  considering  the  fact  that  the  forward  and  back- 
ward transmission  through  turbulent  atmosphere  are  somewhat 
correlated  in  the  case  of  corner-cube  reflector,  whereas  the 
transmission  of  incident  and  reflected  beams  through  turbulent 
atmosphere  will  be  statistically  more  independent. 

With  regards  to  the  peak  frequency  of  the  normalized 
power  spectrum  of  irradiance  fluctuations,  the  experimental 
results  may  be  explained  by  considering  the  fact  that  since 
the  beams  diffracted  by  small  eddies  tend  to  go  outside  of 
the  light  receiver  in  the  folded  path,  the  lower  frequency 
components  due  to  the  larger  eddies  will  be  emphasized  in  the 

3 

folded  path  experiraent  than  in  the  direct  path  experiment. 

In  case  of  corner-cube  reflector,  the  high  frequency  compo- 
nents due  to  the  small  eddies  will  be  scattered  more  than 
single  plane  mirror  reflection  and  its  high  frequency  irra- 
diance fluctuations  will  be  smoothed  out. 

The  frequency  bandwidths  of  irradiance  fluctuations  in 
folded  path  cases  may  be  explained  by  considering  the  Doppler 
shift,  because  the  positively  Doppler  shifted  components  will 


be  further  shifted  to  the  higher  frequency  and  the  negatively 
shifted  components  >rill  be  further  negatively  shifted  in  case 
of  a plane  mirror  reflection.  VIhereas,  in  the  corner-cube  re- 
flector, the  frequency  shift  of  the  reflected  beam  will  be 
in  the  opposite  direction.  (See  Fig.  6 and  7).  In  other 
words,  the  positively  shifted  components  in  one  way  will  be 
negatively  shifted  in  return.  Hence,  the  over-all  frequency 
shift  will  be  smaller  in  corner-cube  reflector  than  the  case 
of  a plane  mirror  reflector. 

By  using  the  angle  of  incidence  fluctuations,  we  should 
be  able  to  formulate  the  power  spectra  of  irradiance  fluctua- 
tions in  corner-cube  and  plane  mirror  reflectors. 
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Statistical  quantities 

^0 

peak 

frequency 

f /f 
o c 

band- 

width 

“hz 

Direct  Collimated  Laser, 
path  length  2.54  m 

0.054 

0.310 

0.886 

0.201 

28 

Comer-Cube  at 

1.27  m 

0.012 

0.278 

0.915 

0.184 

26 

Plane  Mirror  at 

L.27  m 

0.007 

0.187 

1.378 

0.197 

38 

Comer-Cube  at 

2.54  m 

0.025 

0.560 

0.401 

0.126 

38 

’lane  Mirror  at 

2.54  m 

0.103 

0.792 

1.394 



0.140 

50 

TABLE  I.  Comparison  of  measured  quantities  for  corner-cube 

plane  mirror  and  direct  path  where  f measure  peak 

o 

frequency  of  normalized  power  spectrum  and  Tatarski's 

critical  frequency  f = Vx  / / 2ti  AL  . 
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Statistical  quantities 

Yo  ^2 

f /f 
o c 

Direct  collimated  laser, 
path  length  = 274  m 

0.383 

0.124 

0.205 

0.886 

Folded  path  via  corner- 
cube,  one  way  path  = 

274  m 

0.383 

0.592 

0.744 

0.699 

Folded  path  via  plane 
mirror,  one-way  path  = 

274  m 

0.989 

0.695 

0.787 

0.801 

Table  II.  Comparison  of  measured  quantities  for  corner- 
cube  plane  mirror  and  direct  path,  where  f^ 
is  the  Tatarr.ki’s  critical  frequency  V •!•  / 
/Tir  XL  . 
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SCHEMATIC  DIAGRAM  OF  EXPERIMENT 


DATA  ACQUISITION  SYSTEM  FOR  OPTICAL  METHOD 


MIRROR  REFLECTION 


FIGURE  6 
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TRANSMITTER  AND  A RECEIVING  TELESCOPE 
WITH  PHOTOMULTIPLIER 


A CORNER  CUBE  AND  FLAT  MI RROR  REFLECTORS  AND 
A HE-NE  LASER  FOR  DIRECT  TRANSMISSION  EXPERIMENT 


